Current carrying Andreev bound states in a Superconductor-Ferromagnet proximity 

system 
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We study the ground state properties of a ferromagnet-superconductor heterostructure on the 
basis of a quasiclassical theory. We have solved the Eilenberger equations together with Maxwell's 
equation fully self-consistently and found that due to the proximity effect a Fulde-Ferrel-Larkin- 
Ovchinnikov (FFLO) like state is realized in such system. Moreover this state has oscillations 
of the pairing amplitude in either one or two directions, depending on the exchange splitting and 
thickness of the ferromagnet. In particular, using semiclassical arguments (Bohr-Sommerfeld quan- 
tization rule) we show that owing to the presence of the Andreev bound states in the ferromagnet, 
a spontaneous current in the ground state is generated as a hallmark of the FFLO state in the 
direction parallel to the interface. We also discuss the effects of the the elastic disorder and finite 
transparency of the interface on the properties of the FFLO state in the system. 

PACS numbers: 72.25.-b, 74.50.+r, 75.75.+a 



As is well known, ferromagnetism and singlet pair- 
ing superconductivity are competing phenomena. Whilst 
an exchange interaction favors parallel spin aligmcnt, 
Cooper pairs must be in a singlet (spin) state. But 
the two phenomena are less mutually exclusive in artifi- 
cially made ferromagnet-superconductor (FM/SC) het- 
erostructures Q]. In such structures the ferromagnetism 
and the superconductivity can coexist near the FM / SC 
interface owing to the proximity effect 0- In the case 
when the normal metal is not a ferromagnet, the prox- 
imity effect has been studied for a long time 0] • By now 
it is rather well understood in terms of the Andreev re- 
flection processes Q. By contrast, the proximity effect in 
FM/SC systems has became a centre of attention only 
recently. It is not only important from a scientific point 
of view, as it allows for study of the interplay between 
magnetism and superconductivity but also from a 
technological one, as it may find applications in magne- 
toelectronics and quantum computing 6J. 

A number of new phenomena has been revealed in 
FM/SC multilayers. The most interesting examples are: 
non-monotonic behavior of the SC transition tempera- 
ture , oscillations of a pairing am p litu de JM , Q j , lid ] and 
the density of states in the FM [HI [H lia T lj f . param- 
agnetic Meissner effect |15j , proximity induced very long 
range triplet superconductivity in FM \\(\ or genera- 
tion of s pon taneous currents in the ground state of such 
systems These unusual properties, associated with 
Cooper pairs in an exchange field, can be explained in 
terms of an phenomenon first identified by Fulde, Fer- 
rell, Larkin and Ovchinnikov (FFLO) [T3. Originally 
it has been studied in a bulk superconductor with the 
exchange splitting. It turns out that, although of great 
conceptual interest, the bulk FFLO state can be real- 
ized only in a very small region of the parameter space 
near the transition to the normal state [17| and usually is 
destroyed when the exchange splitting E ex is larger than 
V2/2A (Clogston criterion) 18], where A is the SC en- 
ergy gap. Moreover the FFLO state is very sensitive to 



both elastic and spin-orbit scattering 19]. The last two 
effects make the FFLO state very difficult to observe 
experimentally in a bulk samples. The situation is much 
more favorable in FM/SC heterostructures where, due 
to the proximity effect, the Cooper pairs can survive even 
if the exchange field in FM is much larger than SC gap. 

According to our current understanding of the FFLO 
phenomenon in a FM/SC structure, when a Cooper pair 
enters the ferromagnet it acquires a center of mass mo- 
mentum HQ = 2E ex /vF |lflj . where vp is Fermi velocity. 
Usually the SC phase changes linearly with distance x 
from the interface, <f> = Q x x, and this results in a os- 
cillatory behavior of the pairing amplitude in the FM. 
It turns out that under certain conditions a 3D-FFLO 
state, featuring a spatial dependence of the pairing am- 
plitude also along the interface, can be realized [TH l20|. 
It is this latter case that we shall deal with here. 

The purpose of the present paper is to demonstrate 
that in a ferromagnetic layer on a superconducting sub- 
strate, for particular values of the exchange splitting E ex 
and the layer thickness <1f , a FFLO like state is realized 
with the pairing amplitude / varying both perpendicular 
and parallel to the FM/SC interface. It will be shown 
that such a ground state supports a spontaneously gen- 
erated current flowing in opposite directions in the FM 
and the SC regions. The existence of this remarkable 
state was first predicted on the basis of a simple lattice 
model |14| . Here we shall deal with the problem by a less 
model dependent, semiclassical approach and address the 
issue of the observability of the phenomenon in the pres- 
ence of disorder within FM layer and at the FM / SC 
interface. 

The system we consider is sketched in Fig ^ It con- 
sists of thin ferromagnet (FM) of thickness dp deposited 
on a semi-infinite superconductor (SC) and bounded on 
the other side by an insulator. In such an L / FM / SC 
quantum well there will be bound states corresponding 
to the closed quasiparticle trajectories (2l| . Each trajec- 
tory consists of an electron segment, e, which includes 
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FIG. 1: The (Bohr-Sommerfeld) semiclassical trajectories 
for quasiparticles which scatter spectacularly at the I /FM 
interface and by Andreev reflections at the FM/SC inter- 
face. Note that the trajectories have particle-like and hole- 
like segments and they imply an electric current, J p in the 
y-direction. Inset: The energy of the Andreev bound states 
associated with spin up electrons moving in positive (+k y ) 
and negative (— k y ) y-direction as a function of 9 for £f/c(f = 
0.425 and <f> = 0. Evidently, some of these bound states are 
at zero energy. 



FIG. 2: The total density of states p{u>) for £ F /d F = 0.425 
with no current flow (solid line) and in the presence of the 
current (dashed line). For J oc d(j>/dy 7^ the zero energy 
state splits and this lowers the energy of the system. Right 
inset: The difference between total energies of the ABS with 
current flow Egs(4>) an d without the current Egs(0) as a 
function of E ex /A (df/Cs = 1). Left inset: The phase of 
the superconducting order parameter as a function of y. 
The slope d<f)/dy implies a supercurrent J s which is carried 
by Cooper pairs from point a to b (see Fig. 0. 



an Andreev reflection at the FM/SC interface and an 
ordinary reflection at the I /FM interface plus a hole 
segment, h, retracing back the electron trajectory (see 
FigHJ- The Bohr-Sommerfeld quantization rule 



p e (u>)dl 



Ph(ui)dl + 8(f> — 7(0;) = 2nn (1) 



gives the energies of the bound states. The first and sec- 
ond terms represent the total phase accumulated dur- 
ing propagation through the FM region from b to a 
(= Adp(uj + uE ex ) j 'v fCos(9)) , 5<fi is the phase difference 
between points b and a, and j(oj) — arccos(u / A) is the 
Andreev reflection phase shift. An example of the 9- 
dependence of the Andreev bound state (ABS) energies 
for £p I dp = 0.425, where £f = HvF/E ex is the FM co- 
herence length, and (50 = is shown in the inset of Fig. ^ 
Clearly, for any exchange splitting E ex 7^ it is possible 
to find such a 9 that the corresponding ABS is exactly 
at zero energy. If there is a large number of such zero- 
energy ABS for some exchange splittings (more precisely 
for £p /dp ratio) the density of states DOS has a large 
peak at the Fermi energy (ep = 0) . Such a situation turns 
out to be energetically unfavorable. There is a number of 
mechanisms which split this peak and thereby lower the 
energy of the system j23| ■ One of these is a spontaneous 
current which 'Doppler' shifts the quasiparticle energies 
by S = evpA y cos(9) [24|. where A y is a vector potential 
in the y direction. 

An example of such a density of states is depicted in 
Fig. [3 where one can see large peak at zero energy (solid 
line). This corresponds to zero SC phase difference be- 
tween points a and b in Fig. ^ namely no spontaneous 



current. In the inset of Fig. [2 the energy difference be- 
tween state with spontaneous current Eqs(4>) and state 
with no current Eqs(0) is shown. The corresponding 
DOS, with the zero energy ABS split, is plotted with 
the dashed line in Fig. [21 

To summarize, we have shown, using a Bohr- 
Sommerfeld semiclassical argument, that for any ex- 
change splitting, E ex , there are Andreev bound states 
at zero energy and for certain E ex the number of such 
states is so large that it produces huge zero energy peak 
in the density of states. Such a peak in turn is split by a 
spontaneous current and this lowers the total energy of 
the system. This current carrying state can be regarded 
as a realization of the FFLO variation of the pairing 
amplitude in the y-direction. So, one can say that the 
system can be switched between ID and 2D FFLO-\ikc 
states as the exchange field or thickness of the ferromag- 
net is changed. In the following we will show that this 
spontaneous current can be also obtained within a self- 
consistent quasiclassical theory. This approach allows for 
a treatment of disorder and finite transparency of the in- 
terface, and hence provides further useful insights from 
the experimental point of view. 

The quasiclassical matrix Eilenberger equation po| 
reads 



VfVg ff (v f ,r) + 



w CT (r)f 3 + A(r) + — (g <y (vp,r)),g a (vp,r) 

ZT 



where 
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and 



Co a (r) = lo + ioE ex + ievi?A(r) 



(4) 



Here T3 is the Pauli matrix, ui — irT(2n + 1) is the 
Matsubara frequency, a (= ±1) labels the electron spin, 
A(r) is the vector potential and (...) denotes averaging 
over directions of the Fermi velocity vj?. The matrix 
Green's function has to obey the normalization condi- 
tion g^.(vp,r) = 1. The exchange splitting E ex is non- 
zero and constant in the ferromagnet only while A(r) 
is non-zero in the superconductor and is calculated self- 
consistently from 



A(r) = C/7rp T^</ (T (v F ,r)), 



(5) 



where we have assumed that the coupling constant U < 
in the SC and = in the FM. p is the normal state 
DOS and T stands for temperature. 

The spontaneous current has to be determined self- 
consistently together with the Maxwell equation (Am- 
pere's law), which couples the electron current to the 
magnetic field. The total current in the y-direction at 
each point x measured from the interface is given by 



J* ot (x) = 2ienpoTj2{v F g<T(vF,x)), 



(6) 



where e is the electron charge, while the Maxwell equa- 
tion (in the Landau gauge) reads 



d 2 A y (x) 
dx 2 



(7) 



with /iQ being the permeability of free space. 

We have solved the Eilenberger equation numer- 
ically along each quasiparticle 2D trajectory using the 
Riccati parametrization (Schopohl-Maki transformation) 
[2(| together with the self-consistency relations ©-Q). 

The most remarkable feature of the self-consistent so- 
lution is that the iterations of the Eilenberger equations 
frequently converge to a solution with a finite value of 
the current even though there is no external vector po- 
tential. The current flows in one direction over the whole 
ferromagnet and flows back on SC side on the scale of 
SC coherence length £5 (see solid curve in the Fig|3J), 
so the total current is zero, as it should be in the true 
ground state. Such a distribution of the current can be 
associated with a quasiparticle current in FM and with 
a supercurrent on the SC side mediated by Cooper pairs 
(see Fig. [TJ. The fact that the current flows over the 
whole FM is due to the extended nature of the ABS. 
There is also magnetic flux associated with such current 
distribution. Typically the spontaneous magnetic field 
produced by this current is of order of 0.1 B C 2, where 
B c2 is the upper critical field of the bulk SC. 

Within the present self-consistent calculations we were 
also able to study the effect of elastic scattering in the 
FM. The disorder in the (s-wave) SC can be neglected 
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FIG. 3: The total spontaneous current as a function of x 
for a number of mean free path values I. Note that disorder 
introduces oscillations of the current. 



due to Andreson's theorem |27| . The current for a num- 
ber of mean free paths is shown in Fig. [3J As we one can 
see disorder introduces oscillations of the current. The 
spontaneous current is proportional to the DOS at the 
Fermi energy and so the oscillations of the current 
are related to the oscillatory behavior of DOS in the dis- 
ordered sample [l2L Il3l ] . In the clean limit the DOS is 
constant in the whole FM. As is well known, this is a 
property of the Eilenberger equations in the clean limit 
[l3| . Moreover disorder also suppresses the current, as 
expected, since it introduces decoherence of electron-hole 
pairs in the FM. Finally if the mean free path / is shorter 
than the FM thickness dp the current is completely sup- 
pressed. If I < dp the Andreev reflected particles cannot 
reach the I /FM interface, which is a necessary condition 
for the formation of the current carrying ABS, because 
they are scattered on the impurities and the electron-hole 
coherence is lost. In this regard the FFLO variation of 
the pairing amplitude in the y-direction is very sensitive 
to the elastic disorder. However in the x-direction FFLO 
state persists until I < £p [Tol |. even if I < dp. 

To take into account the effect of specular reflections at 
the FM / SC interface we adapt the approach proposed 
by M. Zareyan et al. [3, where a certain probability 
distribution was associated with each semiclassical tra- 
jectory (for details see Ref. 01 )• I n Fig- E]we show the 
current for two different transparencies < r\ < 1. As 
we would expect transparency 77 < 1 suppressed the cur- 
rent because it suppresses Andreev reflection processes 
and at the same time introduces normal (specular) re- 
flections at the FM / SC interface. But even for r\ = 0.3, 
in the present case, we still get the current. It turns out 
that r\ has similar influence on the properties of the sys- 
tem as E ex does. As we can read from the inset of Fig. 
|1] it changes the period of oscillations of the pairing am- 
plitude. Moreover if we changed 77 only, for certain its 
values we get solutions with a current flowing, whilst for 
others solutions with no current. For the above set of pa- 
rameters we have found that current flows in the regions 
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FIG. 4: Effect of finite transparency 77 on the spontaneous 
current E ex /A = 0.2. 77 = 1 means a reflectionless interface, 
while 77 = a perfectly reflective one. Inset: The correspond- 
ing pairing amplitude normalized to its bulk value. Note that 
77 changes the period of oscillations of f(x). 



where 0.7 < 77 < 1 and 0.2 < 77 < 0.35. So we can also 
switch between the ID and 2D FFLO state by changing 
the transparency. 

Finally we note that so far we took no account of the 
demagnetization field B ( i due to the ferromagnet. In the 
limit of a zero film width this is zero. For finite thickness, 
dp , and magnetization, M, (which in the Stoner model 
|28| corresponds to E ex = IM, where I is a fenomeno- 



logical parameter) we estimate Bd to be considerably less 
then B sp due to the spontaneous current. Furthemore it 
should be stressed that in all our calculations the magne- 
tization was constrained to point in the z direction. Thus 
the direction of the spontaneous current was determined 
by the condition that B sp is parallel to M. In a more 
general theory where Andreev orbits also occur in the 
x-z plane and spin orbit coupling is taken into account 
these issues would need to be reexamined. Also the above 
calculations were two dimensional, but preliminary stud- 
ies of a 3D system indicate that there are no qualitative 
changes when orbits in the x-z plane are included. 

In summary we have demonstrated that under certain 
conditions the ground state of a I / FM / SC trilayer fea- 
tures a spontaneous current flowing in opposite directions 
in the FM and SC layers. We argued that this state can 
be viewed as a 2D FFLO proximity state and hence the 
observation of the above current would be a decisive proof 
that the surprising behavior of such heterostructures is 
governed by the FFLO phenomenon. We also showed 
that this state persists only in the clean limit where the 
mean free path is longer than FM thickness and inves- 
tigated the effect of low transparency of the interface on 
the observibility of the ground state current. 
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